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Abstract 



r— H We study ground states of a hybrid system consisting of a polymer and an attractive nanowire by means of computer 
<0 simulations. Depending on structural and energetic properties of the substrate, we find different adsorbed polymer 
^1 conformations, amongst which are spherical droplets attached to the wire and monolayer tubes surrounding it. We 
construct the complete conformational phase diagram and analyze in more detail particularly interesting polymer-tube 



ON 



conformations. 
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1. Introduction 

The study of the interaction between organic and inor- 
ganic matter, or in other words, the behavior of organic- 
inorganic systems, generates fascinating findings with po- 
tential for novel applications in bio- and nanotechnology. 
One of the basic steps in the understanding of such sys- 
tems is the study of the adsorption of soft materials like 
polymers at inorganic matter like solid substrates. In 
the past, numerous computational studies gave general in- 
sights in the adsorption behavior of polymers on planar 
surfaces [H-Q- A particularly surprising fact for example, 
predicted by computer simulations and verified by experi- 
ments recently, is that a single specific mutation in a short 
peptide can substantially change the binding behavior to 
semiconductor substrates @, Q . 

A special class of such hybrid systems are nanotubes 
or nanocylinders interacting with polymers. Carbon nan- 
otubes, for example, are themselves quite interesting nanos- 
jtructures with surprising electronic and mechanical prop- 
erties [8], but nanotube-polymer composites promise to 
enlarge the number of possible novel applications dramat- 
ically, for example in photonics and molecular sensor tech- 
nologies [q|, [Io[- Theoretically, experimentally and com- 
putationally well studied is the wetting of cylindrical sub- 
strates by liquids or polymer droplets. This transition can 
be described by the crossover of barrel-like and clamshell- 



like droplets lH-[l3|. In another study, the adsorption 
behavior of individual polymer chains on nanotubes has 
been studied, where a helical-like winding of flexible and 
semi-flexible chains around the tubes was found [14 1. 
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In contrast, we will here develop a general picture of 
the adsorption behavior of polymers at nanowires depend- 
ing on the properties of the substrate [15]. For this pur- 
pose, we apply a model, where the effective thickness and 
the attraction strength of the linelike substrate are vari- 
able parameters. The above mentioned transitions and 
adsorbed polymer structures are included as special cases 
in this picture. 

2. Model and method 

In our study the polymer is represented by a coarse- 
grained off-lattice bead-stick model, i.e., monomers do not 
have any inner structure and are connected by stiff bonds. 
The polymer is embedded into a three-dimensional simula- 
tion box which includes an attractive thin string pointing 
into the z-direction. The chain is not grafted to the string 
and may move freely. The monomers interact with each 
other via a standard Lennard- Jones potential 



VLj(Vij;e m ,cr m ) = 4e n 
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where is the geometrical distance between two monomers 
i and j and e m and a m are set to 1, such that ^Lj(^min = 
2 1 / 6 ) = — 1. As a remnant of the origin of the model [l6[ 
and in order to facilitate future enhancements and the 
comparison with previous studies, we introduce a weak 
bending stiffness, i.e., the polymer is not flexible in a strict 
way, but may be considered to be flexible in practice: 



^bend(cOS^) = K(l- COs6^) , 



(2) 



where Q{ is the angle defined by the two bonds at monomer 
i and the bending stiffness parameter n is here set to 1/4. 
The interaction between monomers and the string is also 
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based on a simple Lennard- Jones potential, but we neglect, 
as usual [H,[l3j], ^ ne internal structure of the substrate, i.e., 
we assume a homogeneous "charge" distribution along the 
z-axis. We hence simply integrate to get 



V; t ring(r z; ;;ef,cr f ) = a V LJ (Jrl ji ^z 2 ;e f ,aA 
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where r z; i is the distance of the ith monomer perpendic- 
ular to the string and the potential is scaled by setting 
a « 0.528 for convenience 15, I?}- ef and <jf are free pa- 
rameters and can be considered as the string attraction 
strength and the effective "thickness" of the string, which 
is proportional to the equilibrium distance of the string 
potential, respectively. The overall energy of the system 
finally reads 



E 
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VLJ (rij)+^2 Vbend(cOS 0Q + ^ ^string (r Z]i ) • 

(4) 

For estimating the ground-state energies, we apply gen- 
eralized-ensemble Monte Carlo methods 18,[l9[. Confor- 
mational changes are forced by applying a variety of up- 
date moves, including local crankshaft and slithering-snake 
moves and global spherical-cap and translation moves [17j . 



3. Results 

We now discuss low-energy structures of the above de- 
scribed system for different values of the string-potential 
parameters <jf and ef. Based on the simulation of more 
than hundred system parametrizations, we construct the 
full conformational phase-diagram which is shown in Fig.[TJ 




Figure 1: The low-energy conformational phase diagram for poly- 
mers adsorbed at nanostrings. From bottom to top, the string at- 
traction strength ef increases, from left to right, the effective radius 
of the string <7f becomes larger. Different monomer colors or shad- 
ings encode different distances from the string. Monomers near the 
equilibrium distance from the string (colored in red) are defined to 
be in contact with it. 





Figure 2: Visualizations of low-energy conformations with 
N = 30, 100, and 200 monomers in phases (a) Gi, (b) B, 
(c) C, and (d) Ge. 



We identify four major conformational phases, which 
we denote Gi, B, C, and Ge. For small values of both 
<jf and ef, i.e., for weak string attraction and small effec- 
tive radius of the string, we find globular conformations 
with spherical symmetry surrounding the string (phase 
Gi). Increasing the string attraction strength, conforma- 
tions stretch out along the string breaking the spherical 
symmetry and barrel-like conformations with the string 
inside emerge (phase B). In the case of very high string at- 
traction we even find monolayer tubes with each monomer 
being in direct contact with the substrate. Due to the finite 
size of the system, these barrel-like structures break when 
increasing the effective diameter of the string and low- 
energy structures become clamshell-like (phase C), i.e., 
we find adsorbed conformations consisting of a few lay- 
ers which are not wrapping the string completely. Finally, 
decreasing the string attraction at this effective radius, 
conformations become spherical droplets sticked to the 
string (phase Ge). Low-energy conformations from dif- 
ferent regions are visualized exemplar ily in Figs. [T] and [2j 
We convinced ourselves by simulating chains with lengths 
N = 30 and 200 that the general, qualitative structure of 
the conformational phase diagram does not depend on the 
actual length of the polymer, see Fig. [2] for examples. Of 
course, details like the exact positions of transitions lines 
may indeed depend on the actual polymer length. 

To define the different phases and the crossover be- 
tween them, we introduce observables showing a peculiar 
behavior at the transitions. In Gi and B, for example, the 
polymer conformations surrounds the string completely, in 
contrast to the structures in Ge and C. For the localization 
for the transition between Gi-H-Ge and B^C, we hence 
look at the opening angle a [l5[ of the polymer confor- 
mation. The value of this angle shows a jump (low values 
for Gi and B, high values for C and Ge) at the crossover 
between these phases, which is shown in Fig [3l See [l5[ 
for more details on the localization of phase boundaries. 

We finally would like to get at a deeper analysis of some 
structures in phase B. At very high attraction strengths, 
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Figure 3: Opening angle, data points with same cr-value are con- 
nected by lines to guide the eyes. See also Fig.[2fa) and (b) for closed 
conformations (a/ir ~ 0) and[2jc) and (d) for open ones (a/ir > 1). 



rection (see Fig. |4]b). In the adf we clearly see the signals 
from two different regions with different chiralities as well 
as from defects in that structure. A detailed, systematic 
analysis of monolayer structures at high string attraction 
strength and different effective string thicknesses is subject 
of ongoing studies 20[ |21|. 
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low-energy conformations become regular monolayer con- 
formations wrapped around the string, i.e., single- walled 
tubes with an ordered arrangement of monomers form. It 
is noticeable, that there is a competition between differ- 
ent chiral angles, i.e., orientations of the wrapping. This 
behavior is of particular interest as it is in a similar man- 
ner known from carbon nanotubes Figure H] illustrates 
the distribution of the chiral or wrapping angles of such 
a monolayer conformation. Therefore we unzip the struc- 
ture, i.e., we project it onto a plane, and measure the 
angular distribution function (adf) of this unzipped struc- 
ture, whereas we define the chiral or wrapping angle fa of 
the ith monomer as the smallest angle between the vectors 
pointing to its neighbors and the vector in the string di- 




Figure 4: a) Angular distribution function of a monolayer structure 
with (jf = 1.569 (also shown as inset picture in Fig. ^ phase B). 
Inset: Unzipped, planar representation. Different colors (red and 
blue, or light- and dark-gray, respectively) represent main regions 
with different wrapping angles, b) Illustration of the definition of 
the wrapping angle. The arrow points in the direction of the string. 
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